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abstract 

Twisted magnetic flux ropes are ubiquitous in laboratory and astrophysical plasmas, 
and the merging of such flux ropes through magnetic reconnection is an important 
mechanism for restructuring magnetic fields and releasing free magnetic energy. The 
merging-compression scenario is one possible start-up scheme for spherical tokamaks, 
which has been used on the Mega Amp Spherical Tokamak (MAST). Two current- 
carrying plasma rings or flux ropes approach each due to mutual attraction, forming 
a current sheet and subsequently merge through magnetic reconnection into a single 
plasma torus, with substantial plasma heating. Two dimensional resistive and Hall 
MHD simulations of this process are reported. A model of the merging based on helicity- 
conserving relaxation to a minimum energy state is also presented, extending previous 
work to tight-aspect-ratio toroidal geometry. This model leads to a prediction of the 
final state of the merging, in good agreement with simulations and experiment, as well 
as the average temperature rise. A relaxation model of reconnection between two or 
more flux ropes in the solar corona is also described, allowing for different senses of 
twist, and the implications for heating of the solar corona are discussed. 
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1. Introduction 

Magnetic reconnection is a process for re-structuring magnetic fields and rapidly 
converting magnetic energy to thermal energy, with important consequences in many 
laboratory, space and astrophysical plasmas PEE!- Merging of twisted bundles of 
magnetic held lines, known as ” flux-ropes” [HE], through reconnection is widespread, 
and has been investigated in several purpose-built laboratory experiments e.g. PEI El 

El mu EH- 

The MAST spherical tokamak [12] - as well as demonstrating the potential for 
fusion energy generation of the spherical tokamak concept, and providing insight into 
many physical processes relevant to conventional tokamaks - has incidentally provided a 
valuable testbed for reconnection studies. In the merging-compression start-up scheme 
(one of several alternative plasma start-up methods), two flux-ropes with parallel 
toroidal current move together and then merge, creating a single plasma torus with 
closed magnetic flux surfaces [T3]. This provides start-up and current drive without 
using a central solenoid, which could be attractive for future fusion devices, as well as 
allowing an experimental study, with good diagnostics, of reconnection and flux-rope 
merging (EH IHJ in a parameter regime more closely resembling the solar corona and 
other astrophysical plasmas [15] than other experiments. 

The dissipation of stored magnetic energy through reconnection ra - which is 
the primary source of energy release in solar flares - provides a strong candidate for 
resolving the mystery of how solar coronal plasma is heated to temperatures of over 
10 6 K [18, [19]. Such reconnection may occur through merging of twisted magnetic flux 
ropes [SO]. In the flux-tube tectonics scenario [ST], a coronal loop may contain multiple 
twisted threads as its held lines are rooted in several discrete photospheric flux sources. 
Adjacent twisted flux tubes may also be created through photospheric motions with 
multiple vortices within a single flux source. 

In the next section, we summarise recent results from single-fluid MHD and Hall- 
MHD simulations of hux-rope merging in MAST, and discuss the heating of ions and 
electrons through reconnection. In Section |3j we outline how relaxation theory provides 
a useful tool for calculating the final state and the energy dissipated, and present a new 
model of relaxation in a tight-aspect-ratio geometry. Implications for heating of the 
solar corona, mainly based on relaxation theory, are presented in Section [4] 

2. Overview of resistive MHD and Hall-MHD simulations 

MAST [12] is a tight-aspect-ratio toroidal device, with typical major and minor radii R = 
0.95 m, a = 0.60 m, plasma current I p = 400 — 900 kA, toroidal held Bt = 0.40 — 0.58 T 
at R — 0.7 m, and peak electron density and temperature n e0 ~ 3 x 10 19 m -3 and 
T e0 ~ 1 keV. The merging-compression scheme is initiated by the production of two 
toroidal hux-ropes with parallel toroidal currents around the in-vessel P3 poloidal held 
coils (see Figure 1). As the current in the coils is decreased, the attraction of the “like” 
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Figure 1. Cartoon of merging-compression formation in MAST [251 114] . with flux 
ropes shown in purple, forming around the P3 poloidal field coils (left panel) and on 
the point of merging (right panel). 

induced plasma currents in the two flux-ropes causes them to separate from the coils and 
move together, subsequently merging through magnetic reconnection into a single flux- 
rope. Recent experiments using merging-compression have produced spherical tokamak 
plasmas with currents of up to 0.5 MA. The merging is associated with rapid heating 
(presumably as a result of reconnection), indicated by measured electron temperatures 
of up to 1 keV and ion temperatures up to 1.2 keV [121 E2j- 

Recently, fully nonlinear and compressible 2D simulations have been performed of 
flux-rope merging in MAST using both single-fluid MHD and Hall-MHD [H]. The Hall 
terms in Ohm’s Law are considered because, for typical MAST parameters, the ion skin- 
depth is a macroscopic length scale (normalised ion skin depth di = 0.145) and is much 
larger than the width of the Sweet-Parker current sheet predicted by resistive MHD. 
The dimensionless Hall-MHD equations [2] used for the simualations are as follows: 


d t n + V ■ (nvj) = 0 

(1) 

d t (nVi ) + V ■ (mqvi + pi + 77 ) = j x B 

( 2 ) 

E = v e x B - — Vp e + ?/j - r/ffV 2 j 
n 

(3) 

d t B = -V x E 

(4) 


(y-l) 1 [d t p + v, • Vp + 7 pV ■ = rjj 2 + 7] H (Vj) 2 - 77 : Vvi - V ■ q(5) 
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Normalisation is with respect to typical values in MAST start-up plasmas prior to 
merging: density n o = 5 x 10 18 m^ 3 , length-scale L 0 — 1 m, magnetic field B 0 = 0.5 
T, and (equal) ion and electron temperatures T 0 = 10 eV; velocities are normalised by 
the Alfven speed no = B 0 (p, 0 nomi)^ 1 ^ 2 . Here, n is the density, j = V x B the current 
density, v, ; the ion velocity (with electron velocity v e = v* — djj/ri), B the magnetic 
field, p = pi + p e the total (sum of the ion and electron) thermal pressure and E the 
electric field. The ion stress tensor is 7q = — z^(VVj + Vvf), and the heat-flux vector q 
has anisotropic form q = — /cIVyT — T where Vm = b(b • V). The coefficients are 
(all normalised): ion skin-depth, di, parallel resistivity r), parallel ion viscosity v t and 
parallel electron, k\, and perpendicular ion, nf, heat conductivities; here, we set the 
normalised values to be kI = 10” 1 and nf = 10 ‘. The final term on the right-hand 
side of equation ([3]) is hyper-resistive diffusion |2j, used only the Hall-MHD simulations, 
which represents anomalous electron viscosity and sets a dissipation scale for Whistler 
and kinetic Alfven waves, with normalised hyper-resistivity rjjj. Normalised Braginskii 
values of resistivity and parallel viscosity in MAST start-up conditions are r) = 1CT 5 
and = 1CT 3 ; but values of ip v t and rjH are varied in order to study the scaling effects 
of collisions. 

Simulations were performed in: cartesian geometry, with invariance in the 
“toroidal” direction, representing an infinite aspect ratio system; and an axisymmetric 
tight aspect ratio toroidal system, invariant in the toroidal </> direction. The initial 
configuration was chosen to represent the instant when the current rings have detached 
from in-vessel coils, resulting in two localised cylinders/tori of toroidal current. The 
toroidal field within these flux ropes was calculated to ensure local force balance, but a 
non-zero attractive force between the “like” parallel currents caused the rings to move 
towards each other and eventually reconnect. In the tight aspect ratio case, a confining 
vertical field was also incorporated. 

Results from a Hall-MHD simulation in toroidal geometry are shown in Figure 2. 
The key results of these simulations are as follows: 

• The flux-ropes eventually merge into a single flux-rope. For the simulations in 
toroidal geometry, the final state is a single torus with nested flux-surfaces and 
magnetic field profiles resembling those found in MAST plasmas. 

• Oscillations (“sloshing”) in the reconnection rate at low resistivity, resulting from 
magnetic pressure pile-up. 

• A much faster reconnection rate in Hall MHD than in resistive MHD. 

• The qualitative nature of the reconnection in Hall-MHD simulations depends on the 
ratio of the collisional current sheet width - determined by the hyper-resistivity - to 
the ion-sound Larmor radius. At higher collisionality there is a broad current sheet, 
whilst at intermediate values secondary tearing creates small islands of poloidal flux. 
At low collisionality, the outflow separatrices open and fast reconnection is attained. 

• Toroidal simulations in Hall-MHD show the formation of a double-peaked structure 
in the radial density profile which closely resembles experimentally-measured 
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Figure 2. Snapshots showing poloidal field (contours) and toroidal current density 
(colour scale) for Hall MHD simulation of plasma merging in MAST. Toroidal geometry 
was used, with rjn = 10 -8 . 


profiles. 

Finally, we note that preliminary work has been undertaken to predict the evolution of 
ion and electron temperatures, using separate equations for ion and electron temperature 
evolution incorporating, respectively, ion viscous heating and resistive and hyper- 
resistive electron heating [16 j. Plots of the evolution of both temperature profiles 
are shown in Figure 3. Ions are mainly heated by collisional viscous heating in the 
reconnection outflow jets, giving a temperature profile which is double-peaked radially, 
consistent with experiment [25]. The simulations predict ion peak temperatures of 
the order of 1 keV, comparable with the observations of Ono et al ra, but higher 
than reported by Ta.nabe et al [22]; the latter experiments were undertaken at low coil 
currents. The predicted electron temperature profiles are dominated by hyper-resistive 
diffusion which is large in the presence of strong gradients of the current density, and so 
is not cospatial with the current density; this is also suggested by the experimental 
profiles. The experimental observation that the electron temperature is sometimes 
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centrally-peaked, and sometimes hollow, may be linked to the “stochastic” appearance 
and ejection of plasmoids in the simulations. Further work is in progress on simulating 
the evolution of electron and ion temperature profiles. 

3. Relaxation and self-organisation in merging flux ropes in spherical 
tokamaks 

Relaxation theory, first proposed by Taylor to explain the field configurations in RFP 
experiments m, has provided a powerful tool for explaining phenomena in various 
magnetic confinement devices [25] as well as astrophysical plasmas. It is hypothesised 
that a disrupted magnetic field in a highly-conducting plasma relaxes to a state of 
minimum magnetic energy whilst conserving global magnetic hclicity K = J y A • B dV. 
The relaxed state is a linear or constant-/! field: 

V x B = fj, B, (6) 

where /! = /! 0 B • j /B 2 is spatially constant. Note that the relaxed state is a special case 
of a force-free, satisfying j x B = 0. 

It was recently proposed [15] that the merging of flux ropes in MAST could be 
modelled as a hclicity-conserving relaxation to a minimum energy state. The model 
assumes two initial force-free flux ropes separated by a current sheet (a discontinuity 
in the poloidal field); these merge into a lower energy state, consisting of a single flux 
rope within the same volume described by (|6]). In order to permit analytical solutions, 
it is assumed that each initial flux rope has a linear force-free field described by (Jb]). 
The initial field has free energy associated with the current sheet between the two flux 
ropes (this is a delta function layer of reversed current, or negative /r), and hence can 
release energy through relaxation to a fully constant-/! state; it is, however, a force-free 
equilibrium. This represents the field at the moment at which the flux ropes have been 
brought together by the attractive force, but have not yet commenced reconnection. 
Browning et al [15] developed this model assuming infinite aspect-ratio i.e. the flux 
ropes are straight cylinders; a Cartesian coordinate system is used with all quantities 
independent of the axial coordinate z. 

Since MAST, and several other laboratory experiments involving flux rope merging, 
have very low aspect-ratio, it is interesting to extend the relaxation model of [15] to 
finite aspect-ratio. We therefore model relaxation of two toroidal flux ropes within a 
cylindrical ’’can” of rectangular cross-section (hence still allowing analytical solutions). 
The spherical tokamak is assumed to be an axisymmetric configuration contained within 
the volume 0 < Z < /; a < R < h in cylindrical coordinates (R,(p,Z). Lengths are 
normalised with respect to the width of the can, b — a (the minor radius is (b — a)/ 2, and 
the aspect ratio is (b + a)/(b — a). Similarly to the infinite-aspect-ratio model, the initial 
state is given by two constant-/! toroidal flux ropes separated by an annular current 
sheet at the midplane Z = 0.5/, where we take / = 2 (in dimensionless units), giving 
each flux-rope a square boundary similar to the earlier infinite-aspect ratio model eg. 
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Figure 3. Ion (top) and electron (bottom) temperature profiles, in eV, at four 
successive times, for Hall MHD simulation of plasma merging in MAST. The top plot 
also shows in-plane ion velocity vectors. Toroidal geometry was used, with t]h = 1CF 8 . 
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The fields (for both initial and final states) are expressed in terms of a flux function 
-0 as B = (l/i?)(«»<90/<9Z, /ii/j, dip/dR) and equation (jh]) leads to the Grad Shafranov 
equation 


d 2 0 

MR 


1 <90 
RdR + 


<9 2 0 

dZ* 


+ /i 2 0 


0. 


(7) 


The solution closely follows [15], converted from cartesian to cylindrical coordinates. 
The boundary condition is 0 = 0b — constant on the boundary, where 0b (effectively 
a normalisation constant for the fields) must be non-zero. (For 0b = 0, the only non¬ 
trivial solutions to ([7]) are spheromak-like eigenfunction solutions for discrete eigenvalues 
of /i [26]). Following 0 and [26], the solution within a annular shell with height L is 
expressed as a sum of the eigenfunctions 


0 = 0b 


E . . f mi Z 

a nm Rr m (R ) sin I —— 


( 8 ) 


where the radial eigenfunctions are given in terms of Bessel functions R, as 


fmix) = Yi(f3 m b) R(/3 m x) - Yi({3 m x)JiU3mb), (9) 

and /3 m is the mth zero of the function r(R) = Yi(bR)R(aR) — Yi(aR)R(bR). Note that 
equation ([8]) automatically satisfies the required boundary conditions. We normalise so 
that the toroidal flux <3> f = f B z dS — J 0 J b a ^dRdZ = 1 (thereby ensuring conservation 

A 

of flux during the relaxation process described below): hence 0b(/z; a, b) depends on /i. 
Substitution into equation ([7]) gives 

X! [(“Tmn + d 2 )Ran m rm ( R/b ) sin (rmZ/b)] = -/a 2 , (10) 

m,n 

where 


7 2 

/ ran 


= (0 m ) 2 + (mr/L) 2 . (11) 

Multiplying by an eigenfunction r m (R/b)sin(mrZ/L ) and integrating over R and Z, 
exploiting the orthogonality of these functions, then yields: 


(Im.n. 


'7 2 

I mn 


= 0, 


fi 2 4 1 - j m b/a 
T 2 nb 1 - fi 
for m or n even; 


for m, n odd; 


( 12 ) 

(13) 


where 


jm — d\ (/3 m 6)/9i (/3 m o 
The fields from equations 


(14) 


12) are used both to represent the two individual flux 


ropes in the initial state (with L = 1/2 = b — a) and the single relaxed flux-rope (with 
L = l = 2(6 — a)). The final value of //, /if, is determined from the initial value in 
each flux-rope, , by the constraint that both helicity, K, and total toroidal flux are 
conserved (where the dimensionless total toroidal flux initially is 2, because of the two 
flux-ropes). A root-finding process is used to find /i/ so as to conserve K/& 2 . The 
helicity is calculated requiring the toroidal component of the vector potential A^ to 
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Figure 4. Poloidal field lines for initial (left) and final (right) states for the relaxation 
model of merging flux ropes in a tight-aspect-ratio configuration with a =0.1, b = 1.1. 
The initial field (left) has fii = 1.25. The final field (right) has n-f = 0.849. In both 
cases, there is also a toroidal field (Bz = frtp), so that closed poloidal flux contours 
correspond to twisted magnetic flux ropes. 


vanish on the boundary so that the gauge-correction term arising from the flux through 
the torus [28] vanishes. It can be shown that this is related to the total magnetic energy, 
W, by 

K = 2fi 0 W/ia - 2irip b $ t , (15) 


where the dimensionless energy is obtained by integration of the fields, after use of 
Fourier identities for the summation over z and considerable algebra following [26], as 


W = — In 
To 


b 2 

- + /i 2 
a 


: E 

oddm 


2(1 - jmb/a ) 2 


+ R 2 ) + krriQrn + S m ) 


,(16) 


where 


kl = /i 2 - /t, (17) 

P m = l/(cos (k m L) + 1) - 3tan(fc m L/2)/ (k m L) + 1, (18) 

Q m = 2tan(fc m L/2)/ (k m L) - 1, (19) 

Sm = {$n - k m srn(k m L)/L)/(cos(k m L) + 1), (20) 


and '0&(/i, a, b ) is given by the normalisation db = 1 (and /i 0 = 1 in dimensionless units). 

Typical initial and final field states in a tight-aspect ratio geometry are shown in 
Figure 4. The energy released during relaxation is calculated as AW — W t — Wf 
where the initial energy W t = 2W(fii m ,a,b,b — a) and the final energy is Wf = 
4W(nf-,a,b,2(b — a)) (the factors 2 and 4 are to account for, respectively, two initial 
flux ropes each with = 1, and a single final flux rope with <F f = 2. I11 Figure 5, 
the final value of /j and energy release AW are plotted against the initial value, for 
both a tight-aspect ratio configuration similar to MAST and a very large-aspect ratio 
case. Note that the latter graph very closely resembles the infinite-aspect ratio results 
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Final Current [Xf, Energy Release/Length A W' vs Initial Current /ij, with a,b=0.1,l.l 



Final Current fif, Energy Release/Length A W' vs Initial Current fXi, with a,b=50,51 



Figure 5. Final fi = jif (solid curve) and energy released, A W' (dashed curve), as 
a function of the initial value Hi (dimensionless units), for the relaxation model in a 
cylindrical configuration with (top) a = 0.1, b = 1.1 and (bottom) a = 50, b = 51. 

ra- Indeed, it can be shown for an aspect-ratio larger than about 10, // and AW are 
almost identical to the infinite-aspect ratio values. However, at lower aspect-ratios, AW 
significantly decreases; so that the previous infinite-aspect ratio model over-estimates 
the energy release for MAST geometry by a factor of around 2. 

Since the toroidal plasma current is I t = /rdb, Fig. [5] may also be interpreted, 
re-scaling the axes by the appropriate fluxes, as showing the dependence of final 
plasma current on initial current in one flux rope. However, the initial state also 
includes a current sheet, and in order to calculate the total initial plasma current, we 
must add the (negative) net current in this sheet, which may be simply obtained as 
h,cs = -(2/Vo) f b a B r (r,0)dr. 

It can be shown analytically that for small /i*, AW ~ fi }, /i/ ~ /r*. Hence the 
final plasma current depends linearly on the initial current; whilst the energy release 
depends quadratically on /i* and hence also on the initial poloidal field. This agrees 
with experimental findings in MAST and TS-3 [T3j. The dependence of fif on [i n shown 
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in Figure 5 appears as a straight line, but the scaling is in fact stronger than linear, 
which becomes evident for larger /i, : as there is a singularity at the first eigenvalue [26]. 
Such large values of current are not relevant to spherical tokamaks, but may be found 
in solar coronal loops (see Section [5]). 

As shown in Section [2| the magnetic energy released by reconnection may both 
heat the plasma and drive plasma flows - but the latter may also be dissipated to 
provide further heating. We obtain an estimate (strictly an upper bound) on the 
average temperature increase - here simply within single-fluid MHD - by assuming that 
the released magnetic energy is converted fully into thermal energy. Taking typical 
MAST parameters and equating the energy released to the gain in thermal energy, we 
estimate an average temperature increase of 100 eV, which is compatible with what is 
experimentally observed - of course, the relaxation model cannot predict any spatial 
distribution of temperature. 

4. Merging flux ropes in the solar corona 

The simulations of merging flux ropes described in Section [2] may be of some significance 
in understanding reconnection in the solar corona. The latter is, like MAST, a low 
fd plasma, and whilst the Lundquist number in MAST (up to 10 7 ) is much lower 
than the coronal value (around 10 12 ~ 14 ), MAST comes closer to this highly-conducting 
regime than other laboratory reconnection experiments. Coronal reconnection is widely 
modelled using MHD, within which framework the classical ’’Sweet-Parker” model 
predicts reconnection rates which are far too slow to explain coronal heating or 
solar flares, suggesting that physics beyond single-fluid MHD may be required [IS]. 
Furthermore, the Sweet-Parker current sheet width in the solar corona is comparable 
with or smaller than the ion-skin depth, also indicating that Hall physics (at least) must 
be taken into account jT5] . Indeed, analysis of regimes of reconnection suggests that this 
simple comparison of length-scales may under-estimate the importance of incorporating 
Hall physics [29]. The features of reconnection dynamics identified through the Hall 
MHD simulations described in [IT] and in Section [2] above should thus be of some 
relevance to the solar corona. Resistive MHD simulations of merging flux ropes in the 
corona, brought together by attraction of like currents, have been been performed by 
Kondrashov et al ffl; recently, a rather different scenario in which merger is triggered 
by kink instability in one flux rope has also been demonstrated [5T] , 

In order to explain solar coronal heating, it is of primary importance to predict 
the dissipation of stored magnetic energy and, as first proposed by [32], relaxation 
theory thus provides a very useful framework. Browning et al [33] used an approximate 
relaxation model valid for fields close to potential, to predict heating rates in a set of 
adjacent twisted flux tubes. The relaxation model of merging flux tubes [15] may thus 
be adapted to model coronal heating. Since coronal loops have only weak curvature, 
and are not complete tori, the infinite-aspect-ratio model, in which the initial state 
comprises adjacent cylindrical twisted flux-ropes, is appropriate. The model has to be 
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adapted for a loop of finite length l, line-tied to the photosphere at its ends: but it can 
be shown that the gauge-invariant definition of helicity used for an infinite-aspect-ratio 
system (periodic in the axial direction) [15] gives the relative helicity appropriate for a 
field with finite length and non-zero normal field at the boundaries. 

We note that the energy release calculated in [T5J may be approximated by 

AW ~ 0.1 tfr 4 B 2 l//j 0 , (21) 

for a loop of radius r, length l and peak axial field B. Using the relation // ps 2 (f)/l 
where </> is the field-line twist-angle </> (valid near the loop axis) [31], and taking typical 
coronal values B = 0.01 T, l = 10 Mm, r = 1 Mm, we find for rather weakly-twisted 
loops with (j) = Ti an energy release of around 10 19 J. Complete conversion of this energy 
into thermal energy gives a temperature increase of around 10' K for a plasma of density 
10 15 m -3 . This suggests that such flux tube mergers - if occurring sufficiently frequently 
- could make a significant contribution to coronal heating, and might be observable as 
microflares. 

Twisted flux ropes in the corona may arise either through emergence from below the 
surface or by twisting through vortical photospheric motions. Both processes allow for 
the possibility that neighbouring twisted filaments may have twists in opposite directions 
(a situation not possible in the MAST experiment). We therefore consider the relaxation 
of two flux ropes with p = ±/ij, hence having the same fields but with twists in opposing 
directions - this is the counter-hclicity case, whereas the case previously discussed is co- 
hclicity. The total helicity is thus zero, and the relaxed state is therefore potential with 
lif — 0. Hence the energy release is significantly greater than for two flux tubes with 
the same twist, as shown in Figure 6. However, the initial configuration has no current 
sheet since the poloidal field is continuous across the interface z = 0; thus, reconnection 
and merging is less likely in this configuration; if it happens at all - perhaps due to an 
external trigger - it should proceed at a slower rate. On the other hand, evidence of 
such counter-hclicity merging is indeed observed in some solar flares [35], and the large 
energy energy release in such events is consistent with our predictions. 

According to the ’’flux tube tectonics” scenario [23], the solar corona may contain 
many twisted filaments, which might be twisted in mixed directions. In reality, the p 
values of such filaments may cover a range of values, but our analytical relaxation model 
allows us to model only the case of an array of filaments for which the magnitudes of 
/i are all equal (although the signs may vary, allowing for the possibility of filaments 
twisted opposing senses). As an initial state, we assume & n = N x N square array of 
flux-ropes each with p = ip,. First, consider the case of four tubes (N — 2), so there are 
three distinct possibilities for the signs of /if all positive, three positive, one negative; 
two of each sign. (The relative positions of the different orientations is irrelevant to the 
energy release in the relaxation model, which is also unaffected by an overall change 
in sign of twist). As expected, the energy release is greatest in the case when equal 
numbers of positively and negatively twisted filaments are present, so that the total 
helicity vanishes and the minimum energy state is potential, see Figure 7. 
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Energy release in the case of the reversed helicity vs. n- 



Figure 6. Energy release as a function of the initial normalised current (/j,i) for a pair 
of flux tubes with twist in the same direction (red) and opposite directions (blue). 



Figure 7. Normalised energy release as a function of the initial normalised current 
(in) for an array of four flux tubes with twists: all in the same sense; three in one 
sense, one opposite; two in each sense. 


We then investigate the effect of varying the number of initial filaments, n. The 
variation of energy release with initial /i is shown in Figure 8. Here, the magnetic flux 
and the total cross-section are kept constant, so that the width of each initial twisted 
filament is a/N. In all cases, the final relaxed state consists of a single twisted flux tube; 
the energy released increases with the number of flux tubes, due mainly to the increasing 
number of current sheets. The scaling of energy release with number of merging tubes, 
for a fixed /ij is found to be very close to linear AW ~ n 1 ' 07 . This exemplifies relaxation 
as an inverse-cascade from small-scale to large-scale structure. 
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A W(mU|)for different numbers of merging flux tubes 



Figure 8. Normalised energy release as a function of the initial normalised current 
(fii) for different numbers of flux tubes twisted in the same sense, initially in a square 
array within a fixed cross-section: n = 4 (blue), 9 (green) , 16 (red), 25 (cyan), 36 
(purple). 


5. Conclusions 

We have described magnetic reconnection during the merging of magnetic flux ropes 
in the MAST spherical tokamak and in the solar corona - in both cases, reconnection 
leads to the formation of larger-scale magnetic structures and the dissipation of free 
magnetic energy. Simulations of the merging using 2D resistive MHD and Hall MHD 
have been performed, showing in detail how the reconnection rate and dynamics depend 
on the collisionality. The simulations predict the formation of a spherical tokamak with 
a single set of closed flux surfaces; the temporal evolution of the density profile matches 
experimental results. The inclusion of the Hall term in Ohm’s law permits much faster 
reconnection. Whilst the ion skin depth is relatively much smaller, in comparison with 
global length scales, in the solar corona, the Hall term is also likely to be significant 
there, and the simulations are indicative of processes which may be occurring in the 
corona. Simulations show that the magnetic reconnection is associated with significant 
heating of both ions and electrons. Future work will model reconnection in merging 
solar flux ropes, taking account of the Hall term. 

An analytical model based on helicity-conserving relaxation to a minimum-energy 
state allows prediction of the final field configuration and the energy release, which may 
be easily applied to different configurations and parameters. The average temperature 
increase can be calculated if it is assumed that all the dissipated magnetic energy is 
converted to thermal energy. We described how a previous model assuming infinite 
aspect-ratio can be extended to finite aspect-ratio geometry, with particular relevance 
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to spherical tokamaks. The infinite-aspect ratio model provides a good approximate 
model, but the energy release for a given initial ratio of current to field is approximately 
a factor of two lower in the finite aspect-ratio model for MAST parameters. 

The infinite-aspect-ratio relaxation model has been adapted to apply to coronal 
loops. Simple estimates of the temperature increase associated with the merging of 
two weakly-twisted coronal loops suggest this process could contribute significantly to 
coronal heating. In the coronal case, adjacent flux tubes may be twisted in the same 
or in opposite senses - the latter case leads to significantly larger energy release, but 
further simulations are needed in future to explore the conditions in which such flux 
tube merging can occur. A single coronal loop may consist of multiple twisted threads. 
Within a given volume, the energy release increases with the number of twisted threads, 
correlating with the number of current sheets. 
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